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Abstract We describe the data processing pipeline 
developed to reduce the pointing observation data of 
Lunar-based Ultraviolet Telescope (LUT), which be¬ 
longs to the Chang’e-3 mission of the Chinese Lunar 
Exploration Program. The pointing observation pro¬ 
gram of LUT is dedicated to monitor variable objects 
in a near-ultraviolet (245-345 nm) band. LUT works 
in lunar daytime for sufficient power supply, so some 
special data processing strategies have been developed 
for the pipeline. The procedures of the pipeline include 
stray light removing, astrometry, flat fielding employ¬ 
ing superflat technique, source extraction and cosmic 
rays rejection, aperture and PSF photometry, aperture 
correction, and catalogues archiving, etc. It has been 
intensively tested and works smoothly with observation 
data. The photometric accuracy is typically ^0.02 mag 
for LUT 10mag stars (30s exposure), with errors come 
from background noises, residuals of stray light remov¬ 
ing, and fiat fielding related errors. The accuracy de¬ 
grades to be ^0.2 mag for stars of 13.5 mag which is the 
5(7 detection limit of LUT. 

Keywords space vehicles: instruments — techniques: 
image processing — techniques: photometric — tele¬ 
scopes — ultraviolet: stars 

1 Introduction 

Lunar-based Ultraviolet Telescope (LUT) is the first 
robotic astronomical telescope deployed on the moon 
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surface. LUT is pl aced inside a c abin of the lander of 
Chang’e-3 mission (|lp et al.ll2Q14j ) which belongs to the 
Chinese Lunar Exploration Program. After the suc¬ 
cessful landing on the moon of the lander in December 
2013, LUT has been working continuously up to the 
present. It works about 12 terrestrial days per month 
in lunar daytime for sufficient power supply. LUT ob¬ 
servation has two key programs: pointing observation 
program and survey program. The pointing program 
monitors brightness of variable targets which are pro¬ 
posed by world-wide astronomy society, including cat¬ 
aclysmic variables, RR Lyrae stars, eclipsing binaries, 
active chromosphere stars, flaring M dwarfs, etc. From 
the first light on Dec. 15, 2013 to Feb. 2015, LUT has 
monitored variable stars for ^800 hours, and surveyed 
an area of ^1600 deg^ around the moon’s north pole. 

Typically, one target is continuously monitored for 
^50 hours, which generates about two thousands of im¬ 
ages and requires about 4.5 Gbytes storage space. Data 
processing is not performed real-time because data are 
obtained after LUT finished its lunar daytime work. 
Data processing of pointing observations should be fin¬ 
ished in a handful of days so as to leave enough data 
processing time for other observation programs. There¬ 
fore, an automatic data processing pipeline is manda¬ 
tory to address the issue of massive data quantity. The 
most parts of the pipeline follow the astronomy data 
reduction routines, while some parts are developed for 
special features of LUT observations. For example, the 
LUT images suffer from significant pollution of stray 
light from scattered sunlight in the cabin and the tele¬ 
scope. Another pollution source comes from the cosmic 
rays in space, which is much more significant than that 
on earch. Furthermore, flat field images taken with 
internal LEDs are not perfect in terms of large-scale 
uniformity. To solve this problem, the large-scale il¬ 
lumination structure, i.e. superflat, obtained through 
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dithering observations, is coupled with the illumination 
corrected and normalized LED flat field. 

This paper describes the data processing pipeline re¬ 
ducing the pointing observation data of LUT. A brief 
description of LUT’s instruments and pointing obser¬ 
vations are presented in Section [2j Details of the pro¬ 
cedures of the pipeline, and the building of LUT flat 
field are described in Section [3l Section 0] describes the 
aperture correction method. The precision of the pho¬ 
tometry after the stray light removing and flat fielding 
are shown in Section [5l 


2 Instrument, Observation Calibration 


LUT is a 150 mm, F/3.75 Ritchey-Chretien telescope 
working at a Nasmyth focus. A flat mirror is mounted 
on a two-dimensional gimbal in front of t he telescope 


apert ure for pointing and tracking (see I Wang et^ 


201 Fig. 1). The mirror can rotate from -28° to 
+13° in azimuth (axis of telescope as zero), and +20° 
to +38° in altitude (horizontal direction as zero). A 
UV-enhanced back-illuminated AIMO CCD is mounted 
on the fucal plane, and UV coating is applied on one 
lense of the field corrector as the UV filter. The result¬ 
ing passband of LUT is about 245-345 nm, peaking at 
250 nm. The CCD pixel scale is 4.76” pixeU^, so the 
exposure area of CCD with 1024x1024 pixels gives a 
field of view (FOV) 1.35°x 1.35°. Two pairs of LEDs 
are installed crosswise (one as backup) on the front in¬ 
side wall of the CCD camera, which can be used to 
illuminate the CCD through a ring-like diffusing glass 
for flat field calibration. The LEDs emit at 286 nm 
and the spectral widths are ^12 nm. For further de¬ 
tails of the scientific objectives, instrumentation, sys- 
te m performance, a nd c alibrations of LUT, please refer 


to 


Cao et al. ( 2011 1 and Wang et al. I (l2ni,^ . 


The pointing observation strategy is described as fol¬ 
lows. A target is placed near the center of the CCD and 
monitored for several observational runs. Each run lasts 
for about 30 minutes, consists of several exposures, and 
has a fixed telescope pointing with respect to the moon. 
During a run, the total shift of the target due to the ro¬ 
tation of the moon is within a region of ^ 50 x 100 pixel. 
The shift of stars in image during each exposure is small 
(within 1 pixel) comparing with their profile widths of 
about 2pixels (see Sect. 13.4p . The next run re-direct 
the flat mirror pointing to make the target return to 
the center. Such a strategy favors stray light removing 
and flat fielding (see Sect. Ej). 

Calibration observations include dark field acquisi¬ 
tion, internal flat field exposures and superflat obser¬ 
vations. LUT can obtain internal flat field images to 


correct pixel-to-pixel nonuniformity in CCD sensitivity, 
making use of its LED lamps. However, the LED illu¬ 
mination is not ideal in terms of large-scale uniformity. 
To correct the large-scale nonuniformity, superflat im¬ 
ages are created employing dithering observation tech¬ 
nique. Before dithering observation, a positional grid 
was designed to sample the large-scale nonuniformity 
structure in the FOV. Usually, a grid size of 7x7 was 
adopted (see Fig.[T]). At each nodal point, the standard 
star was observed for about 20 times. 

Superflats are not created making use of the sky 
background, because the atmosphere on moon is ex¬ 
tremely tenuous so there is no sky background available, 
and also because the ecliptic light is not in LUT’s avail¬ 
able sky area. But after all, the stay light suffered by 
LUT would certainly contaminate the flat fields. There¬ 
fore, the superflat is actually uncovered using standard 
stars as uniform light sources. The dithering obser¬ 
vation was carried out at most once for each month, 
depending on both the Chang’e-3 and LUT operation 
plan arrangement. They had been carried out in Jan¬ 
uary, June, August and December 2014, and in January 
and May 2015. Flat field correction in each month made 
use of the superflat of the adjacent month. 



Fig. 1 A typical dithering observation sampling grid of 
the flat field taken in June 17, 2014. The target star was 
HD 152303 with signal-to-noise ratio (SNR) 66.8 measured 
in 2xFWHM aperture radius. XCENTER and YCENTER 
are column- and row-direction coordinates of the star on 
every image frame. 
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3 Pipeline Description 


The data processing pipeline is developed for LUT 
pointing observations data reduction, focusing on ob¬ 
taining catalogues and light curves of targets of inter¬ 
est. The pipeline consists of several procedures, includ¬ 
ing overscan correction, stray light removing, astrom¬ 
etry, flat fielding, source extraction, cosmic ray rejec¬ 
tion, aperture radii determination, photometry, aper¬ 
ture correction and catalogues archiving procedures, 
which are summarized in Table [H Each procedure is 
described in the following subsections. The pipeline 
is developed with SExtractor, IRAE and PyRAE 
(IRAE wrapped in Python) programs, and Python 
packages including NumPy, SciPy, AstroPy, Py- 
EITS and astLib. 

LUTs raw data are originally obtained from the Data 
Management Subsystem (DMS) of Ground Research 
and Applicatio n System (GRA S) of GE-3, in Level OB, 
binary format (|Tan et al. 2Q14h . After data delivering, 
all data are converted to PITS (Elexible Image Trans¬ 
port System) format. The EITS headers record the 
LUT observation modes and instrument working status. 
The keywords in the headers include “TASKCODE” - 
the type of observation task, “IMAGETYP” - the type 
of images, “ELE” and “AZIMUTH” - the pointing co¬ 
ordinates of the flat mirror (for details see Table [2]). 


3.1 Stray Light Removing 


The LUT detection suffers from stray light problem 
caused by sunlight being scattered by the cabin and the 
telescope. The strength and pattern of the stray light 
are varying, with ADU counts from a few thousands in 
common cases to a few tens of thousands, depending on 
the angular distance between the flat mirror pointing 
and the sun. Eortunately, in most cases the variation 
of stray light evolves very little in subsequent images 
taken within ^0.5 hour, so a method had been devel¬ 
oped to remove the stray light from those images, as is 
described below. 

Eirstly, all images are preprocessed through overscan 
correction and trimming to size of 1024 x 1024. Then, 
images are grouped according to their head keywords 
“TASKCODE”, “ELE”, “AZIMUTH”, and the expo¬ 
sure time of each group should be within less than 1900 
seconds. Eor each given image, a specific stray light 
pattern is derived from the other images in the same 
group except itself through image combination using 
the IRAE “median” algorithm. According to the point¬ 
ing observation strategy, the stars positions on succes¬ 
sive images always have slight shifts of a few pixels, so 
during the “median” combination all celestial sources 


are rejected and a stray light pattern is hence left in 
the combined image. Also contained in the combined 
image are the underlying bias level and dark current 
counts. Then, each image subtracts the combined im¬ 
age, thereby removing the stray light pattern, bias, and 
dark counts. A group commonly contains 15-30 images, 
so the combination can give high SNR stray light tem¬ 
plates. Thus, the pattern removing procedure can be 
considered not to induce extra noise to images. Our 
stray light removing procedure also removes underly¬ 
ing bias and dark current in the mean time. Eigure [2] 
gives an example of stray light pattern and the result 
after its removal. 


3.2 Astrometry 


The astrometry is performed by cross matching the star 
distrib utions on LUT images with the Tycho-2 cata¬ 
logue ( Hog et al. I2QQQ[) . which have been trimmed to 
match LUT’s total available sky area. Out of every 
frame 5-10 bright stars are extracted using SExtrac¬ 
tor and fast photometry of the stars are performed. 
The geometrical distribution and the measured LUT 
magnitudes of the bright stars are then cross matched 
to the Tycho-2 catalogue. The star coordinates in the 
catalogue have been transformed to the current epoch 
and have also been corrected for precession. Tycho-2 
optical magnitudes have to be converted to LUT NUV 
magnitudes for brightness matching. They are firstly 
transformed to standard Vega B, V magnitudes through 
the transformation relationshi ps provided by the Hip- 
parcos and Tycho Catalogues (|ESAI [l99^ . Then, the 
Vega B, V magnitudes are used to calculate theoretical 
LUT AB m agnitudes through t h e stel lar atmosphere 
model from Castelli and Knruczl ( 2004 ) for a series of 
different spectral types (details are described in a paper 
of Han et al. in preparation). Eor each matched star, 
in addition to the cross matching radius constraint, the 
difference between the measured and calculated LUT 
magnitude is used as further constraint and is required 
to be less than ^2 mag. If the cross matching successes, 
the world coordinate system (WGS) (J2000) as a result 
is written into the EITS header. If the matching fails, 
the pipeline gives up the current frame and jumps to the 
next one. If more than 20% of total frames fail in the 
matching, the automatic data processing pipeline work¬ 
ing with WGS terminates and an alternative pipeline 
working without WGS will be carried out. The accu¬ 
racy of astrometry is typically about 1”. 


3.3 Elat Eielding 


The flat fielding procedure makes use of both the in¬ 
ternal flat field from internal LEDs and superflat from 
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Table 1 Outline of Data Processing Pipeline 

Data Processing Outline 

0. Data Preparation 

1. Overscan Correction & Image Trimming 

Stray Light Removing 

2. Image grouping according to “AZIMUTH”, “ELE” and time period 

3. Combine images of a group adopting “median” algorithm to make stray light pattern template 

4. Each image subtract its stray light template 

Calibration 

5. Astrometry 

6. Elat Eielding With Rectified Elat Image 

7. Source Extraction with DETECT_THRESH=2 and DETECT_MINAREA=4 
Profile measurement and Clipping 

8. Clip objects in margins of images and clip ELONGATION>2 objects 

9. Measure Moffat profile EWHMs for every objects in an image 

10. Keep objects that have 1.3<MEWHM<3.2 

11. Determine typical MEWHM for an image group and assign to EWHMmed 

Photometry 

12. Aperture photometry with apertures radii in units of EWHMmed 

13. PSE photometry 

14. Aperture Correction 

15. Calculate center and corners’ J2000 coordinates and write into EITS header 

16. Data and Catalogues Archiving, light curve output 


Table 2 FITS header keywords recording the LUT observation modes and instrument work status. 


Keyword 

Value 

Comment 

TASKCODE 

“Initial” 

Default observation task 

“Pointing” 

The pointing observation task 

“Astrometry” 

Astrometry calibration observation task for the telescope 

“Survey” 

The survey observation task 

IMAGETYP 

“Object” 

Image of celestial objects observation 

“Zero” 

Image of zero calibration 

“Dark” 

Image of dark calibration 

“Elat” 

Image of flat field calibration 

ELE 

number in arcsec 

Altitude angle of the flat mirror 

AZIMUTH 

number in arcsec 

Azimuth angle of the flat mirror 



Fig. 2 The effect of stray light removing. Left: An example of original LUT CCD image in size of 1072 x 1027; middle: 
The derived stray light pattern of the left image made from images of its host group in size of 1024 x 1024; right: Image 
after stray light removing in size of 1024 x 1024. 
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dithering observations (see Sect. [2j). The internal flat 
field has been processed to retain only the CCD pixel- 
to-pixel response nonuniformity, filtering out the low- 
frequence structure, which produces illumination cor¬ 
rected and normalized LED flat field. The dithering ob¬ 
servation of a single standard star produces a positional 
sampling grid of LUT’s large-scale response nonunifor¬ 
mity. These nonuniformity structures are recorded in 
both the background and the stars in each image. Af¬ 
ter stray light removing, the information of large-scale 
response in background is removed, but is left in stars’ 
fluxes. Since each standard star is used as an “invari¬ 
able” light source, its fluxes at different positions in 
the FOV should exhibit the structure of large-scale re¬ 
sponse nonuniformity. The flux counts at these grid 
positions are measured through aperture photometry, 
before which, flat fielding using the internal flat field im¬ 
age is performed. Then, the fluxes of the grid is fitted 
by a two-dimensional, second-order polynomial func¬ 
tion and the superflat is created (see Fig. [3] top). The 
final fiat field that is used in the pipeline is the product 
of image multiplication of the superfiat and the illu¬ 
mination corrected and normalized LED flat field (see 
Fig. [3] bottom). 

3.4 Source Extraction, Brightness Profile 
Measurement and Cosmic Ray Rejection 

Sources in each frame are extracted using S Extrac¬ 
tor by certain criteria to conveniently obtain their 
CCD X-Y coordinates. The key extraction criteria 
are summarized in Table [3l The selected parame¬ 
ters for SExtractor output are XJMAGE (object’s 
barycenter position along X-axis), YJMAGE, ELON¬ 
GATION (shape parameter, major-axis/minor-axis), 
BAGKGROUND, etc. 

In order to determine the aperture size for pho¬ 
tometry (see Sect. 13.5|) . brightness profile measurement 
and cosmic ray rejection are performed for each im¬ 
age. First of all, sources in edge regions are cleaned 
from the extracted catalogue because they may be sig¬ 
nals arose just by instruments. Then, FWHMs of 
sources profiles are measured by the “psfmeasure” task 
in IRAF.noao.obsutil package. X- and Y-axis posi¬ 
tions in the output file from SExtractor are used 
as inputs of the task, and Moffat profile widths (de¬ 
noted as MFWHMs) are required as output. In LUT 
images, cosmic ray events are comparable to celestial 
sources in number. The event rate of cosmic rays 
is estimated to be ^4 cosmic rays per second in the 
FOV, corresponding to ^2 cosmic rays per second per 
square degree. We adopt the cosmic ray rejection cri¬ 
terion as MFWHM<1.3. The criterion is determined 


through sources identification in some typical cases, in 
which cosmic rays are identified by correlation of suc¬ 
cessive images. Figure 0] shows the MFWHM distri¬ 
butions of the identified stars (red) and cosmic rays 
(blue) in a typical case. Although a few cosmic rays 
are blended with stars in terms of their profile, the dis¬ 
tributions of the two populations are obviously seper- 
ated, which enable us to reject most cosmic rays in 
terms of MFWHM^l.3. Objects with MFWHM>3.2 
are deemed as extended sources or cosmic ray clumps, 
so they are also clipped off from the extracted cata¬ 
logue. Residual cosmic rays are further rejected with 
the shape criterion of ELONGATION>2, if the target 
of pointing observation is not a binary or blended stars. 
The clipping criteria are summarized in Table [3l 



Fig. 4 An example of histogram statistics of MFWHMs of 
celestial objects (red) and cosmic rays (blue) for a single im¬ 
age. A clipping line is set at MFWHM=1.3 to remove most 
cosmic rays from the extracted catalogue. A ID gaussian 
fit is performed (black line) to the histogramic distribution 
profile after 1.3<MFWHM<3.2 filtering. 


3.5 Aperture and PSF Photometry 


Both aperture and PSF photometries are performed 
for all the extracted objects after the cosmic ray re¬ 
j ection . The AB magnitude system of Oke and GnnnI 
(Il983h i s adopted following S loan Digital Sky Survey 
(SDSS, iFnkugita et al. 19961). Galaxy Evolution Ex¬ 
plorer (GALEX. iMorrissev et al. 2007 ). etc. LUT mag¬ 
nitude is defined as 


R^lut = R^o,LUT — 2.5 log /lut (1) 

where mo, lut is the zero point magnitude of LUT, /lut 
is the flux density in ergss“^ cm“^ Hz“^). Zero point 
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X (pixel) 



Fig. 3 Flat field creation. Top: two-dimensional, second-order polynomial fitting to superflat; bottom left: an image of 
the two-dimensional, second-order polynomial fitting to superflat (a); bottom middle: the processed internal flat field image 
only retaining pixel-to-pixel response nonuniformity (b); bottom right: final flat field image for calibration = (a)x(b). 


Table 3 Source extraction criteria for S Extractor and clipping 


Source extraction criteria for SExtractor 

Criteria 

SExtractor parameters set up 

Detection threshold relative to background RMS 
Minimum number of connected pixels above threshold 
Background mesh size 

DETECT.THRESH = 2 
DETECT.MINAREA = 4 
BACK SIZE = 64 

Clipping Criteria For Output Parameters of SExtractor And IRAF.psfmeasure 

Criteria 

Parameters set up 

Clip edge region of CCD columns 

Clip edge region of CCD rows 

Clip stretched objects 

Clip cosmic rays and extended sources 

20 < XJMAGE < 1020 

5 < Y.IMAGE < 1020 
ELONGATION < 2 

1.3 < MFWHM < 3.2 
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magnitude of photometry is calibrated for LUT system 
as 


^o,LUT = 17.52 zb 0.05. 


( 2 ) 


Details of LUT p hotometry calibration are given by 
Wang et al.l (|2015 b The corresponding signal-to-noise 
ratios (SNRs) are calculated as: 


SNR^. = 


iFr, , , X 


(3) 


(IRAF.apphot.phot HELP document) where F^. is the 
total number of counts excluding background in aper¬ 
ture n, G is the gain of CCD (electrons per ADU), 
A is the area in aperture ri in square pixels, a is the 
standard deviation of the background which mainly in¬ 
cludes readout noise, stray light noise, bias and dark 
counts noise, local flat fielding noise, etc., A^sky is the 
pixels number of background. 

The aperture photometry is performed with a series 
of aperture sizes, whose radii are lx, 1.5x, 2x, 2.5x, 
3x, and dxFWHMmed, where FWHMmed is denoted as 
the stars’ typical FWHM. The FWHMmed is calculated 
for each image group which has formed in the stray 
light removing procedure. The calculation method 
of FWHMmed is US follows. Firstly, for each image, 
MFWHMs that satisfy 1.3<MFWHM<3.2 are plotted 
in a histogram with bar width of 0.2 pixel. Secondly, a 
typical FWHM for each image is obtained by fitting the 
MFWHM distribution with a one-dimensional Gaus¬ 
sian function and deriving the Gaussian peak value 
(denoted as MFWHMpeak)- Instead of a PSF mea¬ 
surement based on the brightest star, this statistical 
method is adopted because there are some cosmic ray 
events whose profiles perfectly mimic that of a bright 
star, which can not be removed through the rejection 
criteria described above. Thirdly, the median value of 
MFWHMpeak of a group of images (15-30 in number 
varies for different pointing observation tasks) is de¬ 
rived and denoted as FWHMmed, and is used as the 
unit of aperture size. The average value (“median” 
algorithm) is adopted here to deal with possible star 
brightness profile variation, although such variation is 
negligible (within ^0.04 pixel through an observation 
task) for most cases. The variation may only be sig¬ 
nificant (in a factor of ^2) in the first observational 
day of each month, if the telescope has not reached 
its designed thermal equilibrium state. Background 
annulus for aperture photometry is set to be 6x and 
8 X FWHMmed as inner and outer radii, respectively. An 
input file of sources positions for each image is used as 


IRAF image cursor, so aperture photometry runs auto¬ 
matically. The photometric error is typically ^0.01 mag 
for LUT band 10 mag stars (30s exposure), which comes 
from image background noises. 

Before PSF photometry, a further clipping based on 
FWHM criterion is carried out to select candidate PSF 
stars. We firstly clipped the extracted sources with 
FWHM<1.4 and FWHM>2.5 pixels. After the clip¬ 
ping, the brightness profiles of the 10 brightest ob¬ 
jects are fitted through minimization. A test work 
of comparing various profile functions, such as Gaus¬ 
sian, Lorentzian and Pennian functions, has been car¬ 
ried out and indicates that an elliptical Moffat func¬ 
tion with a fixed /3 parameter of 1.5 provides best fits 
more frequently. After the PSF model establishment, 
PSF photometry is carried out within circle radii of 
3 X FWHMmed- The inner radius of the annulus used to 
determine background level is 5x FWHMmed, and the 
annulus width is 2 x FWHMmed- The effects of PSF fit¬ 
ting for single and double star are shown in Figure [5l 


4 Aperture Correction 


The correction of aperture effect is applied to the mag¬ 
nitudes measured by the aperture and PSF photometry 
to obtain the magnitude measured in an “infinite” size 
aperture, which is considered to embrace the total flux 
of a source. The aperture correction is performed by 


7T?/'j^-^cor — TThf^ | 


(4) 


where is the aperture correction factor for aper¬ 

ture ri, and m^.^cor are magnitudes before and 
after aperture correction, respectively. The values 
of the correction factors are determined through the 
“curve of growth” method. The curve of growth is 
derived by performing aperture photometry for stan¬ 
dard star HD 185395 with a dense sampled aperture 
radii series. The 22 sampled aperture radii range from 
0.5 X FWHMmed to 25 X FWHMmed- The standard star 
has 43 frames of exposure, and 43 curves of growth are 
obtained through aperture photometry. Mean values 
of magnitudes at each aperture radii are calculated to 
derive the mean curve of growth. Figure [6] {left) shows 
the mean curve of growth in terms of the magnitude 
offsets relative to magnitude measured within radius of 
14xFWHMmed- The solid line in Figure [6] (/e/i^) is the 
modeling approach to the mean curve of growth, which 
can be expressed by the following function: 


Am{r) = < 


1 — 2.5l0g(l — e 1.93-2.29r+l.4r^-0.27r^ )^ 

(r^2) 

1.22 - O.lr + O.OlGr^ - 0.00086r^ 

(r>2) 











Fig. 5 PSF fitting for single and double star. The hot colored surfaces illustrate source fluxes, and the blue wired frames 
illustrate PSF models to the sources. Left: PSF fitting with Moffat function for a star with SNR~31; right: PSF fitting 
with Moffat function for a star with SNR~43 in a binary star BW Dra. 


where r is the desired aperture radius to which the pho¬ 
tometry is corrected in units of FWHMmed- 

Aperture correction factors are calculated as mean 
differences between magnitudes in every apertures and 
the magnitude in the aperture from which magnitude 
differences converge to zero and also converge to their 
errors. The criteria of converge are 

^ 0 ( 5 ) 

I ^ — mrj (^) 

where is the magnitude measured in aperture n of 
the mean curve of growth, 

is the error of —nir^. Figure [6] (n^M) 
shows the magnitude differences and their errors ver¬ 
sus aperture radii. From SxFWHMmed on and up 
to IbxFWHMmed aperture radius, the magnitude dif¬ 
ferences converge. A mean value of magnitudes in 
10 X, 12 X, and 14xFWHMmed aperture radii is calcu¬ 
lated and used as the magnitude of total flux. Mag¬ 
nitude offsets between lx, 1.5x, 2x, 2.5x, 3x, and 
4xFWHMmed aperture photometry magnitudes and 
the total-flux-magnitude are used as aperture correc¬ 
tion factors, whose values are listed in Table IH 

5 Pipeline Performance 

We have carried out some test work to assess the ac¬ 
curacy of the data processing pipeline. Since the flat 
field images are made from the processed internal fiat 
field images retaining pixel-to-pixel nonuniformity, and 
the superflat reflecting large-scale nonuniformity, they 
would not contain medium-scale (tens to one hundred 
pixels) structures, which may bring flat fielding related 
errors to photometry. The influence of medium-scale 


Table 4 Aperture correct factors (in magnitude) and their 
errors. 

Ti ArUr^ 


(FWHM^ed) 


1.0 

-0.348 


7.3E-4 

1.5 

-0.150 


6.5E-4 

2.0 

-0.086 


6.5E-4 

2.5 

-0.059 


7.2E-4 

3.0 

-0.043 


7.3E-4 

4.0 

-0.024 


7.3E-4 


structures can be tested by high frequency positional 
sampling observations whose targets go across the im¬ 
age frame. At each position, aperture photometry is 
performed and magnitude is obtained to find out the 
medium-scale structures effect. Figure [3 shows the re¬ 
sult of such test carried out in June, observing stars 
HD 204770, HD 205022 and HD 203711, where XCEN- 
TER means X-axis positions of the targets, MAGlut 
means LUT magnitudes measured in 3xEWHMmed 
aperture radius. Standard deviations of magnitudes are 
0.020 mag for HD 204770, 0.022 mag for HD 205022 and 
0.021 mag for HD 203711. Such dispersion is mostly 
caused by medium-scale nonuniformity, which is esti¬ 
mated to be ^0.02 mag after deducting errors of image 
background noises (^0.01 mag for a 10mag star). The 
error increases to ^0.2 mag for stars of 13.5 mag with 
30 s exposure, and this corresponds to the 5cr detection 
limit of LUT. 
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Fig. 6 Left: LUT point source photometric curve of growth. Magnitude of 14xFWHMmed aperture is adopted as 
the reference line, and aperture axis maximum range for clear show; right: |mri^i,med — '^^ri,med| {green eireles) and the 
corresponding errors {red squares). 



Fig. 7 Photometry test for integrated effect of stray light 
removing, ffat fielding and aperture photometry. 


6 Summary 

We describe the data processing pipeline reducing the 
pointing observation data of LUT. The pipeline per¬ 
forms stray light removing, astrometry, flat fielding em¬ 
ploying superfiat technique, source extraction, source 
profile measurement, cosmic ray rejection, aperture and 
PSF photometry, aperture correction, catalogue archiv¬ 
ing, and outputs light curves. The pipeline has been in¬ 
tensively tested and works smoothly with observation 
data. The photometric accuracy is typically ^0.02 mag 
for LUT 10mag stars (30s exposure), with errors from 
all that of background noise, residuals of stray light re¬ 
moving, and fiat fielding. The accuracy degrades to be 
^0.2 mag for stars of 13.5 mag which is the 5cr detection 
limit of LUT. 
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